Can we detect Hot or Cold spots in the CMB with Minkowski Functionals? 
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In this paper, we investigate the utility of Minkowski Functionals as a probe of cold/hot disk- like 
structures in the CMB. In order to construct an accurate estimator, we resolve a long-standing 
issue with the use of Minkowski Functionals as probes of the CMB sky - namely that of systematic 
differences ("residuals") when numerical and analytical MF are compared. We show that such 
residuals are in fact by-products of binning, and not caused by pixelation or masking as originally 
thought. We then derive a map-independent estimator that encodes the effects of binning, applicable 
to beyond our present work. Using this residual-free estimator, we show that small disk-like effects 
(as claimed by Vielva et al. [1, 2]) can be detected only when a large sample of such maps are 
averaged over. In other words, our estimator is noise-dominated for small disk sizes at WMAP 
resolution. To confirm our suspicion, we apply our estimator to the WMAP7 data to obtain a null 
result. 
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I. INTRODUCTION 

Cosmological inflation has become a promising paradigm to solve many problems of standard cosmology in a unified 
framework [3] , providing a compelling model to explain large scale flatness and homogeneity of the universe in addition 
to the existence of small density perturbations which form the seeds of structure. 

Among the possible inflationary scenario is false vacuum driven eternal inflation [4]. An inflaton potential with 
multiple minima allows for the nucleation of bubbles with smaller vacuum energy density. This might occur due to 
tunneling via the Coleman-DeLuccia mechanism [5] or via the collision of such bubbles [6, 7]. Pressure differences 
between the bubble walls force the bubble to expand, quickly approaching the speed of light. In this picture our 
observable universe can be thought of as a bubble residing among a multitude of bubbles in a 'multiverse'. 

However, bubbles inevitably collide and recent work shows that such collisions may leave observable imprints in 
the CMB. A generic prediction is that a past collision on our bubble universe by another bubble will leave a cold or 
hot circular disk - regions where the mean temperature is statistically different - on the CMB sky [8-12]. In addition 
to such shift in mean temperatures, the CMB may exhibit additional polarization modes [13] in such regions and 
perhaps lead to anisotropic large scale galaxy flows [14]. There have been some claims in the literature regarding 
the existence of such a spot in the CMB, the so-called "cold spot" [15-17] detected using wavelet analysis, although 
such claims have been challenged [18, 19]. More recently, a model independent pipeline was constructed to search for 
such signals using causal boundaries and found several possible hints of such features [20, 21]. In this paper we will 
attempt to search for such a signal using a different statistic - Minkowski Functionals. 

Since spots in an otherwise smoothed gaussian sky are topological in nature, this suggests the use of statistical 
descriptors which are well suited to quantify morphological properties of the temperature fluctuations. Minkowski 
Functionals are exactly such tools - they are morphological statistics on smooth maps. While they are widely used 
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in image processing (e.g. [22]) as such, they were first used by cosmologists to look from deviations from Gaussianity 
of the perturbations in large scale structure [23-25] and the CMB [26-34] . 

In this paper we apply MFs to the search for disk-like structures in the CMB - structures expected if our present 
"bubble universe" has had the (mis)-fortune of colliding with another bubble in the distant past. In section II we 
introduce MFs for scalar fields on the sphere. We review the standard results for Gaussian random fluctuations and 
show how these results change in the presence of a bubble collision. Section III begins with a brief introduction to 
the numerical method we have set up to simulate Gaussian maps - in particular we derive analytical formulas to 
remove numerical "residuals" introduced by binning and derive a general map-independent residual-free estimator. 
In Section IV, we propose our analytical model for disks in the CMB and apply our estimator to constrain both 
simulated collisions maps and the WMAP7 data. In the last section we summarize and conclude. 



II. MINKOWSKI FUNCTIONALS OF SCALAR FIELDS ON THE SPHERE 



A. General definition 



MFs characterize the morphological properties of convex, compact sets in an n dimensional space. A property is 
considered to be morphological when it is invariant under rigid motions, i.e. translations and rotations. Hadwinger's 
Theorem under some weak assumptions ensures that any morphological property can be expanded as a linear com- 
bination of n + 1 MF. On S 2 there are 2 + 1 = 3 MFs which, up to normalization, represent the area, circumference 
and integrated geodesic curvature of an excursion set. For a given threshold is, the excursion set and the boundary of 
the excursion set of a smooth scalar field u on the sphere are defined by 

Q u := jx e § 2 | u{x) > t/j , (la) 
dQ v := Ix G § 2 | u(x) = . (lb) 

The first MF [44] vo (v) is the area fraction of Q u , given by 

— ( dtt&{u-v) , (2) 
4tt J S 2 



where is the Heaviside step function. The second MF is proportional to the total boundary length of Q v 



V ^ V )~TT l dl = TT I dQS(u-u)\Vu\ . (3) 

167r JdQ„ 167r JS 2 

Here 5 is the Delta distribution and |Vu| is the norm of the gradient of u. Finally, the third MF is the integral of the 
geodesic curvature k along the boundary 

«a(f) := 3-2 / 3-2 / dSl5(u-v)\Vu\K, (4) 

The geodesic curvature, k, describes the deviation of the curve 7 from being geodetic. For a unit speed curve, i.e. 
I7I = 1, it is defined through 

k := |Vyy| , (5) 

where represents the covariant derivative along the tangent vector 7 of the curve. Thus k vanishes if and only 
if 7 is geodetic. However, for the numerical calculation of v^, it is convenient to express k in terms of u. To do so, 
one can use the fact that u does not change along 7 and thus du (7) = which implies that 7^ = e^V v u. Upon 
normalization, this can be used in equation (5) to yield n in terms of the metric and derivatives of u. 



3 



B. MFs for a Gaussian random field and a superimposed collision 



The expectation values of the integrals in eqs. (2) - (4) for a Gaussian random field m g are given by [35, 36] 



tf(«0 :=<,o G H) = Ml-erf(^-=|) 



*f» := <«?(„)) = ^cxp , (6b) 

where erf is the Gaussian error function erf(x) := ^= dtexp (—t 2 ). The mean, variance and its moment are given 
by 

(i:=(u G ), *:=(u 2 g )-h 2 , t := \ (|V UG | 2 ) . (7) 

The variance a and the variance of the gradient can also be expressed through the Gaussian angular power spectrum 
Cf by 

- = i$> + i)cf> r=±±V + l)C? l «±}±. (8) 

However, note CMB power spectra are often truncated at large I - hence it is preferable to compute these quantities 
directly from the maps using eqn. (7). 

Assuming that the primordial spectrum is completely Gaussian and the power spectrum is isotropic, we pose the 
question: If a "cold" or "hot" spot exists in the CMB (whatever the origin) - how well can we distinguish such a spot 
from the complete gaussian sky with MF? By "hot" or "cold" spot, we mean a circular region of size A in the CMB 
with a uniform temperature shift ST over the actual mean temperature /i G = M — A/(Att)ST and the usual power 
spectrum of anisotropics, where ^ is the average temperature of the CMB, and we have assumed a sharp cut-off at 
the boundary. Furthermore, the actual variance of the unaffected region of the sky o~f is related to the disk properties 
by oq = a — Aj (47r) (1 — A/ (An)) ST 2 , where a is the variance of the whole sky calculated assuming that no such 
disk exist. In general, the "hot" or "cold" spot can be fairly complicated in structure, with gentler boundaries and 
non-uniform temperature shift. In this work, we will ignore such complications and consider a single sharp boundary 
region with uniform temperature shift. Such regions are predicted by studies of cosmological bubble collisions [12]. 

In addition, we would like to point out some subtleties when using MF to test for the presence of disk-like structures. 
First, MF are very weak probes of the two-point correlation function, so one in principle can have a completely gaussian 
map, yet is somehow correlated pixel by pixel - imagine for example, an omnipotent hand rearranging all the cold 
pixels in a gaussian map into a disk, then we will not pick up such a feature. Since we are assuming that the underlying 
power spectrum is isotropic, such magical rearranging does not occur. This illustrates one of the advantage of using 
MF over "local" statistics - we will not mistake fortuitous (but random) correlations as a true feature. 

Second, imposing a gaussian disk with different mean on the sky, effectively renders our map to become bi- 
distributional, i.e. a histogram of the pixels will reveal two different gaussian distributions with different means 
but identical variances, zero skewness and kurtosis. Hence the sky becomes non-Gaussian. Nevertheless, this is a 
very specific form of anisotropic non-Gaussianity, which cannot be described by a regular higher-point correlation 
function such as the bispectrum or trispectrum [34] which is isotropic by construction. Thus, in principle we can hope 
to distinguish such anisotropic non-Gaussianities from those generated during primordial inflation [37-40]. We will 
show later that they possess a very distinctive signature. 

For MFs of a Gaussian map containing a disk, we propose the decomposition 

f A \ A 
Vi(v,fi,a,T) := ( 1 - — J vf(v,^ G ,o- G ,T G ) + —vf (v - ST, /i G ,cr G ,r G ) + dA,(v, ST, shape) , i e {0, 1,2} . (9) 

The first two terms are area weighted MFs of pure Gaussian fields. The first term represents the part of the sky 
which is unaffected by the collision. The second term corresponds to the MFs of a Gaussian field in A with the mean 
temperature shifted by ST. The third term dAi stands for the boundary effects of the transition region where the 
temperature drops from ST 0. Note that, within this decomposition, information about the shape of the collision 
region is entirely contained in dAi since the first two terms solely depend on the constant temperature shift ST and 
the area of the collision region A. 
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III. NUMERICAL MINKOWSKI FUNCTIONALS, RESIDUALS AND COSMIC VARIANCE 
A. Expected Residuals in the numerical calculation of MFs 

For the numerical extraction of MF for a collision map, we employ the HEALPix suite of tools [41]. We use 
HEALPix to generate a full sky map of temperature anisotropics from a given power set of a; m 's using the power 
spectrum derived from the five year WMAP data [42] . Our MF code is a straightforward numerical calculation of the 
integrals in cqs. (2)-(4) and closely follows the prescription in [27] as follows. 

Given a pixelated map with field value u-iXi), we can calculate its first and second partial derivatives at each pixel 
- in HEALPix this is done in (I, m) spherical harmonic space. The numerical MFs Vi are computed via a sum over 
all pixels of the respective integrands 

2o (v>%j) ■= &(u — v) , (10a) 



l t (u, Xj) := i<5 (u-v)y/ u 2 ^ + v? v , (10b) 
X 2 ( V| Xj ) := -6 (u v) -.-^ ^ , (10c) 



where 



u.$ := o®u , u ;v := - — -cLu , u ;w := , ? - o u + — — -o®u , 

1 „ _ cos^ „ d 



(11) 



smv sm -Q ov 



Summing over all pixels we obtain the numerical MFs 

i_ N < 



V i {u):=—-^2l i (u,x j ). (12) 



The integrands X\ and I2 involve the delta function which is numerically approximated through a discretization of 
threshold space in bins of width Av by the stepfunction 

S N (x) := (Av)" 1 [9 (x + A1//2) - 6 (x - Ai//2)] . (13) 

This numerical prescription was proposed for use in the analysis of CMB maps by Schmalzing and Gorski [27], 
and then implemented in subsequent analysis by various authors [30-34]. In these works, it was noticed that even 
when given a completely Gaussian map, the MF calculated from the above prescription possess systematic "residuals" 
which is resolution and sample size independent. This is usually attributed to pixelation effects, masking effects and 
other intangibles. The standard procedure is to estimate the residuals by generating a large sample of Gaussian maps, 
calculating the MF numerically, average the MF over the large sample, and then subtracting from it the analytic 
prediction Vi(y), 

Ai(v):={Vi(y))-Vi(u). (14) 

This residual Ai(y) is then subtracted from all other numerically calculated MF, the implicit assumption being that 
it is the same even when the underlying map is non-Gaussian. However, as we will show below, this residual is in fact 
a byproduct of binning - it is the substitution of a delta function to the discrete delta function eqn. (13). Further, 
we will show that the residual is map dependent and that it scales as bin-size squared (Ai>) 2 in the leading order. 
In the following, we derive an exact expression for the residuals applicable for a general underlying map, and show 
numerically that the calculated residuals arc reproduced. 

First, notice that the procedure corresponds to the replacement of the delta function S(u — v) in the integrals 

/oo 
du 5{u - v)vi(u) , ie{l,2}, (15) 
-OO 

with the numerical delta function 5n(u — v). However, this results in 

V i {v)=v i {v)+R? v {v), (16) 
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FIG. 1: Average numerical MFs for a Gaussian random field generated from the power spectrum derived from 
the five year WMAP data at -ZV s id = 512 with -d s = 1° smoothing (red dots), compared to the expectation value vf 

(blue line) as given in eqs. (6a)- (6c). 
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FIG. 2: The figure shows the difference Af (eqn. (19), red) of the numerical and the analytical MFs and the leading 
order residual R^q with Av/^/a ~ 0.46 (blue) for an average over 256 realizations at A s id = 512 and i? s = 1° 
normalized to max \ v!f\. The standard deviation of these maps is ^fa ~ 0.065mK. 



with residuals defined through 



v+Av/2 



if-Au/2 



du Vi{u) - Vi(v) . 



(17) 



As it does not depend on the actual functional form of Vi, this is a general result and represents a generic map- 
independent (but binsize dependent) MF estimator. 

In the special case of a Gaussian random field, Vi = v^, we expect that the numerical computation of the MFs 
yields, 



^ gM= 8A^V2 



erf 



v — n+ Ai>/2 



'2a 



erf 



v - [i - Avf2 



i 



(2tt) 3 / 2 a Av 



exp 



{v-ii-Avjiy 
2^ 



12a 
exp 



[y - (i + Av/2Y 
2^ 



(18a) 
(18b) 



We tested our numerical implementation of MF against a completely Gaussian map, for which the integrals (2), (3) 
and (4) are known analytically. Figure FIG. 1 shows a comparison of the numerical MFs with the expectation 
values (6a), (6b) and (6c). The average of the difference 



ArH :=V?{y,n,a,T)-ti?{v,n,o,T) , i € {1,2} , 



(19) 
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FIG. 3: The figure shows the difference A^ (eqn. (19), red) of the numerical and the analytical MFs and the leading 

order residual (blue) for a single realization taken at 7V s id G = 512 smoothed to § s — 1° at the binwidths 
Av/y/a ~ (0.92, 0.62, 0.46) normalized to max|t>|*|. The upper (lower) panel is Af (AJf). As shown here, smaller 
binsize leads to smaller residuals but increased noise - for a single realization we find that a binsize of Av j ^fa ~ 1 

works well for 7V S id = 512. 



between numerically extracted MFs and the respective expectation value, as shown in Figure FIG. 2, is in very 
good agreement with i?fg when a sufficiently large number of realizations is considered. All maps are corrected for 
the numerical fluctuation in the mean such that /j, = O (10~ 18 ). The difference Af is normalized to the maximal 
amplitude in the MFs 

max\v$\ = l, max|^| = iy^, max | v° | = ± I . (20) 

As we emphasized previously, the calculation of the residual can be done for any underlying smooth map. For 
completeness and future reference, we will derive the residuals for hierarchically non- Gaussian maps in an Appendix. 

An expansion of eqns. (18a) and (18b) around Av — shows that the leading order terms are proportional to 
(Av) 2 /a. This fact may suggest that a smaller binsize is always better. However, smaller binsize means that each 
bin contains fewer pixels and hence an increase in the inherent noise per bin. In figure FIG. 3 we show that, for a 
single realization, a binsize of Av j ^fa = 0.9 is a good compromise for iV s id c = 512. Another way of beating down 
the noise is to increase the number of pixels, either by increasing the resolution of the map, or average over a large 
sample of maps. Figure FIG. 4 shows a comparison of the residuals R^q with the difference Af at the prospective 
Planck resolution A s ;d c = 2048 for a single sample at different binsizes. With the number of pixels per bin increased 
a smaller binsize can be used without adding noise. 

However, as we would want to apply our prescription to actual data, we ultimately want to be able to extract 
accurate MF from a single map. As we shall see, the consideration of only one realization will turn out to be a 
difficult stumbling block (due to cosmic variance) in our attempt to constrain disk-like structures in the sky. 
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FIG. 4: The figure shows the difference Af (eqn. (19), red) of the numerical and the analytical MFs and the leading 

order residual R^q (blue) for a single Gaussian realization taken at -/V s id e = 2048 without smoothing at the 
binwidths Av/yfa ~ (0.60, 0.46, 0.31) normalized to max \vf\- The upper (lower) panel is Af (Aff ). The map has 
variance a = (O.lllmK) 2 . With increased number of pixels that comes with higher resolution, we can use smaller 

binsizes while keeping the pixel noise tolerances manageable. 

B. Remaining difference in MFs of Gaussian maps 

The residual effects that originate in the numerical implementation of the delta function have been analyzed in 
detail in the last subsection. Henceforth, we will use a suitable binwidth in the calculation of MFS for any given map 
to obtain well converged residuals. Their subtraction then allows for an efficient removal of the effects of the discrete 
delta function. Consequently we are interested in further effects that may cause a difference between the numerical 
MFs and their analytical expectation. Therefore we investigate the difference 

AfM^W-^M+iZ^M] , <& = 0, (21) 

that remains after the residuals have been removed. Examples of this remaining difference in the MFs of one realization 
of the Gaussian field are shown in figure FIG. 5, while the average of the differences are shown in figure FIG. 6. 
The difference in the MFs of a particular realization differs from sample to sample and therefore has a random 
character. Moreover, for most samples the shape of the difference Af appears to be dominated by ± {\fad v ) vf 
times some random prefactor of the order (O (0.1)) 3 . The averages of the difference, (A^V converge to a curve that 

is approximately given by (O (0.1) V^,) vf and are thus much smaller than the random fluctuation in the MFs of 
a single sample. We point out that this difference in a single sample depends on the resolution (though only weakly 

2. /3 

oc iV sido ) and does not depend on the binwidth Av. 
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FIG. 5: The figure shows the difference Af, eqn. (21), of the numerical and the analytical MFs with residuals 
subtracted for three different realizations taken at iV s ido = 512 without smoothing at the binwidth Av/^/a = 1 and 
normalized to max \ vf I . Though the differences Af are different in each sample they appear to be dominated by 
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FIG. 6: The figure shows the difference Af, eqn. (21), of the numerical and the analytical MFs with residuals 
subtracted, averaged over 1000 realizations taken at N^c = 512 without smoothing at the binwidth Av/^/a = 1 
and normalized to max \vf\. The averages (Ap) appear to converge to a curve that is approximately given by 

(O (0.1) \fadv) vf. However, when compared with the upper figure, it turns out that the fluctuations for a single 
realization are about an order of magnitude larger than the average. 



IV. ANALYSIS OF COLLISION MAPS 



A. Expected disk signal in MFs 

In this section, we describe our method of constructing a sample gaussian map with a superimposed disk. 
Let A = 2ir (1 — cos$£>) be the the area of the disk of angular size with temperature difference ST. We linearly 
add this into a map of gaussian temperature anisotropics uq by 

u D {■&) = ST ■ 9 (#£> - tf) . (22) 

In practice this is done in spherical harmonic {aim) space via the sum a\ m = af m + a^ ml where af m is the gaussian 
spectrum and 



af m = ST J ^-j (Pi-. 1 (cos0 D ) - P l+1 (cos# D ))S m0 , (23) 

is the disk spectrum. Both spectra are cut off at the some high l max > 1000, which has to be high enough to ensure 
that the steepness of the step function is preserved [45]. While in principle, there exist a small contribution from 
the boundary region of the disk, as the transition region is very small due to the steepness of the step, the signal 
associated with the gradient dAi is highly suppressed and hence we neglect it from now on. Thus we will henceforth 
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FIG. 7: Average numerical MFs of a Gaussian random field with a superimposed disk of §d = 60°, ST = S^/ac a t 
N s ide = 512 (red dots), compared to (blue line) as given in cq. (24). 



consider 

Vi(u,iJ,,a,T) := (l - — j vf{v,i± G ,a G ,T G ) + — (i/ - ST, fi G ,a G ,T G ) , (24) 

as a sufficiently accurate approximation to the expectation values of MFs from a collision map. However, the downside 
is that this also implies that we cannot access information about the shape of the boundary because it is contained 
in the boundary terms [46]. Numerically extracted MFs for a disk of temperature difference ST = 'iy/ac and opening 
angle -do = 60° are shown in figure FIG. 7. 

Notice that equation (24) is invariant under the simultaneous replacement of ST —> ST = —ST and A — >• A = 4ir — A. 
This is a simple reflection of the fact that a "hot" spot of tempcratue ST and size A in a Gaussian field with expected 
temperature [i G may equally well be regarded as a "cold" spot of temperature —ST and size 4ir — A within a Gaussian 
field of mean temperature n G + ST. This degeneracy can be circumvented by restricting the consideration to disk 
sizes A < 2tt. 



B. Relative amplitude of the disk signal 



In this subsection we investigate which kind of disks one can hope to detect using MF, by comparing their form 
expected from a Gaussian field with those expected in the presence of a disk with temperature difference ST and 
opening angle $d- Consider the difference 

Avi(v, (i, a, t, ST, A) := v { {v, fi, a, r) - vf{v, fi, a, r) , (25) 

where \i, a represent mean and variance of the field and r is the variance in the gradient field of a given temperature 
map which is supposed to contain a disk. The disk parameters are constrained by a G = a— A/(4-7r)(l— A/ (An))5T 2 > 0. 
In figure FIG. 8 we have shown the differences Avi for a disk with ST j %fo = 0.7 and do = 50°. Unfortunately, the 
shape of the difference is very similar to d^vf for a wide range range in the {ST, $£>}-parameter space. This means 
that the difference Af that occurs in the numerical MFs of a Gaussian field can potentially be mistaken as the 
presence of a disk. Only when the temperature difference clearly obeys ST/ yfa > 1 and the disk covers a significant 
fraction of the sky a distinctive signature manifests itself in MFs, hence the temperature difference and size of a 
disk must be large to be detectable with MFs. Moreover, the intrinsic variation in a single outcome of numerically 
generated Gaussian map at WMAP resolution is large enough to mimic the presence of a prominent hot or cold spot 
in the MFs of the map. In other words, MFs are not a very sensitive tool when it comes to the detection of disks in 
the CMB. 

This issue is elucidated in figures FIG. 9 and FIG. 10 where we show the normalized L 2 -norm of the difference 

Ai(y) :=Vi(i/)- [Vi(v, n, a, t)+R? 1 '{v, h,*,t)] , (26) 
R-"(v,H,o-,t) := (\ - ^\ Rf (v, m ,a G ,T G ) + ^-Rf {v - ST, ^ g ,o g ,t g ) , Rq U = , (27) 
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FIG. 8: The Figure shows the normalized difference Avi for ST '/ ' \fo = 0.7 and &d = 50° with the n,a and r taken 
from the realizations presented in figure FIG 5. Notice that Avi is very similar in shape as d^vf so that the 
difference that occurs in Af may potentially be mistaken as the presence of a disk of this amplitude. 
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FIG. 9: The figure shows lnLf for a single Gaussian realization and its minimum (red dot) at WMAP resolution 
(left, iV s ide = 512, without smoothing). The ST = line is degenerate in space. 
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FIG. 10: The figure shows lnLf for an average over 1000 Gaussian realizations and its minimum (red dot) at 
-^Vgide = 512 without smoothing. The ST = line is degenerate in $£> space. 



i.e. 

L\ := """^<ff'» 2 , ( 28) 
max \ Vi\ 

for single Gaussian realizations at WMAP and for an average over 1000. Figure FIG. 11 shows the £ 2 -norm for 
prospective PLANCK resolution A^ide = 2048. The minimum is indicated by the green dot and therefore means a 
best fit of eq. (26) to the data. The apparent presence of a disk in the case of single realizations the fluctuations in 
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FIG. 11: The figure shows mL 2 for a single realization and its minimum (red dot) at N s ia c = 2048 without 

smoothing. The ST = line is degenerate in $d space. 



the MFs of the Gaussian field, as shown in figure FIG. 5, is due to the fact that their shape is very similar to the 
shape in the difference Avi, figure FIG. 8 so that eqn. (24) allows for a good fit for the data. Figure FIG. 10 shows 
that upon averaging over a larger number of samples the fluctuations in the MFs decrease, cf. figure FIG. 6, and the 
best fit essentially resembles the null result. Notice the degeneracy of the ST — line in ■d^ space. 

In figures FIG. 12 and FIG. 13 we show an example of fake collision data for ST/yfa = 1 and i9u = 50°. As in 
the Gaussian case, the remaining difference in the MFs has a severe effect on Aj so that the minimum of its L 2 -norm 
(red dot) is not to be found at the input values (green dot). However, when we average the MFs taken from many 
realizations, the remaining difference in the MFs decreases and the minimum of the L 2 -norm of Aj is very close to 
the actual input parameters. We conclude that we cannot detect disks with ST/^/a < 1, even if they cover a large 
fraction of the sky. Only large disks with $£> = O (10°) with temperature difference ST / y/a > 2 for which the main 
contribution to the MFs lies clearly outside of the Gaussian, cf. figure FIG. 7, can be detected with certainty. The 
main drawback to the use of this MFs algorithm is the remaining difference Aj which results in a bad signal to 
noise ratio for these disks. As this difference depends only weakly on the resolution we do not expect a significant 
improvement from PLANCK data. 



C. Application to WMAP seven year data 

In the recent literature on the CMB cold spot, see e.g. [15-17], it has been argued whether the occurrence of such a 
spot is a likely feature of a Gaussian random field and therefore is regarded "generic" , or whether it is a non-Gaussian 
feature the origin of which is related to a thus far unknown physical mechanism. 

In this context we ask what can be inferred about the presence of spots from the MF statistics. Therefore we have 
used the same map making procedure as [43], Viclva et al. [1] and Zhang and Hutercr [18] in their analysis of the 
CMB cold spot. We compute the temperature u of the fiducial map at Xi by the sum 

u( Xi ) = £^>f . (29) 

of temperatures of each individual differential assembly r € {Q1,Q2,V1,V2,W1,W2,W3,W4} weighted by w r (i) = 
N r (i)/a^, where N r (i) is the number of effective observations at pixel number i and ay is the noise dispersion of the 
respective receiver. The maps are added at a resolution of N s ^ e = 512 and smoothed with <& s = 1°. The a; m 's are 
extracted before the KQ75 mask is applied. The MFs are then computed by summing only over the unmasked pixels, 
i.e. 

with Wj = 1 when the pixel is not hidden behind the mask and otherwise. 

It is clear that MF does not have the sensitivity to pick up the small signal as seen by [1], i.e. a ST = — 0.016mK 
at 5 degrees at roughly 3er, since in this regime the signal is smaller than remaining noise of a single realization 
as described in section (IV B). Indeed, fitting the co-added map into our estimator, we obtained a "best-fit" with 
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FIG. 12: The figure shows In if for a single realization with ST = yfa and do = 50° (green dot) at iV s ido = 512 
without smoothing. The corresponding best fit value is indicated by a red dot. 
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FIG. 13: The figure shows lni? for an average over 2048 realizations with 6T = yfa and i?tj = 50° (green dot) at 
^Vgidc = 512 without smoothing. The corresponding best fit value is indicated by a red dot. 



temperature difference ST ~ — 0.063mK and opening angle "&d — 35°, which is clearly a fit to noise and hence is not 
physical. 



V. SUMMARY & CONCLUSIONS 



Motivated by recent work on cosmic bubble collisions and their potentially observable signatures in the CMB, we 
studied the utility of MFs for their detection. 

In order to do this, we resolved the long-standing issue with the MF "residuals" - systematic differences between 
analytically and numerically computed MF which arc independent of map resolution and sample sizes. We show that 
these residuals are in fact a result of finite bin-sizes, and not caused by pixclation, masking or other intangible effects 
as originally suspected. We derive a map-independent analytical formula to characterize these residuals at all orders, 
allowing one to convolve effects of bin-size into the MF estimators. 

After removal of these residuals, we find that the remaining discrepancies between the analytic estimates and the 
numerical MF are of order O (lO -3 ). This discrepancies is proportional to the number of pixels of the map and the 
number of sample sizes, indicating that we are approaching the limit expected from random noise alone. Unfortunately, 
as we demonstrated in the text, this noise has a characteristic that is roughly similar to the expected disk signal, and 
hence severely limits our ability to probe small disk signals. 

We apply our residual-free MF estimator to the investigation of Gaussian temperature fluctuations containing a 
superimposed collision signal. To characterize the signal-to-noise of our estimators, we generated collision maps by 
modeling the signal as a uniform shift of mean temperature within a circular spot (a disk) in an otherwise Gaussian 
field. We find that our least-squares fitting procedure accurately reproduces the underlying signal only when a large 
number of realizations of maps are averaged over. For a single WMAP and PLANCK resolution map we are able to 
recover the result only when a highly prominent disk, with \ST\ > 2^/ctg and > 40° is present. This is unfortunate, 
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as it means that MF are intrinsically too noisy to be able to distinguish cold and hot spots in the CMB for small sizes 
as suggested by [1]. In order to confirm our suspicion, we apply our prescription to WMAP7 map and find that we 
do not recover the latter's conclusions. 
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Appendix: Residual removal in hierarchically non-Gaussian maps 

In section III we mentioned that in [30-32, 34] the residual that is removed is an average over the residuals of 
Gaussian maps, cf. equation (9) in [32]. However, the total residual from the delta function (for i = 1,2) also contains 
non-Gaussian terms which we like to present here. We refer to the notation used in [34] and note that their definition 
involves the normalzed threshold Vm which relates to the threshold used in the present work as vm = (y — )i)/^/g, 
while their variance gm translates as gm = yfa and g\ = s/2t. Henceforth we use x = vm and gm = c. Accordingly, 
the MFs can conveniently be expressed as 

V i (x)=A i exp(-x 2 /2)v i (x), ie{0,l,2}, (31) 

where Ai are constants that are related to the two point correlation function. In the hierarchically non-Gaussian case 
the functions Vi(x) can be expanded in powers of the variance g as 



(32) 



The individual terms are given by (x) = Hi-i(x) and 



= ^H i+2 (x) - ^H t (x) i{i -^ SlI Hi _ 2 (x) , (33a) 

v W( x ) = ^-H 5 (x) + ^H 3 (x), (33b) 

vf\x) = gif 6 (z) + ^=^W) 1 [k, + |tf ) H 2 (x) , (33c) 

vi 2) (x) = ^H 7 (x) + K ~ 2 f Sl H 5 (x) - i (kj + ifiTfi-n) H 3 (x) \ (k u + ^iSn) H x (x) . (33d) 

The Hi{x) are Hermite polynomials defined by H n (x) = exp (x 2 /2) {—d/dx) n exp (— z 2 /2), where for n = —1 we have 

H-\{x) = exp (x 2 /2) y/ir/2 (l — erf (x/V2)) and the constants S and K are the skewness and kurtosis parameters 
and their ri-th derivatives S n and K n . 

The residual that is to be expected in the second and third MF from the numerical approximation of the delta 
function via a stepfunction of width Av, equation (13), is given by equations (16) and (17). That is, the numerical 
MFs yield 

Vr m (x) = V i (x)+R^{x), ie{l,2}, (34) 

which has the residual contribution 



R^(x) :- 



1 

Ax 



x+Ax/2 

dxVi(x) 

x—Ax/2 



Vi(x) . (35) 
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This motivates the introduction of the new estimator Vi(x) which corrects for the residuals by 

Vi{x) := Vr^ix) - R? "(i) , ie{l,2}. 

In analogy to eqn. (31) 

a px+Ax/2 

RfrJx) := / dx exp (-i 2 /2) »r } (5) , with Ax := Av/a . 

Ax J x -Ax/2 

Consequently, the respective residuals are given by 



[exp(-i 2 /2)iV 2 (x)]==B£Z(x), 



Ax 

A 

Ax 



s 



iSi 



Ax 

M 

Ax 



exp {-x'/2) ( -H i+1 (x) - -fHi-xlx) - 
( ^H 5 (x) 



K-S Si 
24~ 



i(i - l)5n 



1 



K - 2S S\ 
24 
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1 



Hi- 3 {x) 
3 



F 3 (5)--UTi + -^ F x (5) 



Km 



H- X {x) 



Hi(x) -q[Ki + -SS u ) H 2 (x) --[Ku + -SiS n 



1 



(36) 
(37) 

(38a) 
(38b) 

(38c) 

1 

(38d) 



where [/ (x)] ± = f (x — Ax/2) — f (x + Ax/2). The residuals of hierarchically non-Gaussian maps up to order two 
arc therefore 



(39) 



n=0 
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